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o 

£N| ' Abstract. In this paper we consider parametric ideals and introduce a notion of compre- 

hensive involutive system. This notion plays the same role in theory of involutive bases as 
the notion of comprehensive Grobner system in theory of Grobner bases. Given a parametric 
ideal, the space of parameters is decomposed into a finite set of cells. Each cell yields the 
corresponding involutive basis of the ideal for the values of parameters in that cell. Using 
the Gerdt-Blinkov algorithm described in [6] for computing involutive bases and also the 
Montes DisPGB algorithm for computing comprehensive Grobner systems [13], we present 
an algorithm for construction of comprehensive involutive systems. The proposed algorithm 
■ has been implemented in Maple, and we provide an illustrative example showing the step- 

' by-step construction of comprehensive involutive system by our algorithm. 



1 Introduction 



, One of the most important algorithmic objects in computational algebraic geometry is Grobner 

basis. The notion of Grobner basis was introduced and an algorithm for its construction was 
designed in 1965 by Buchberger in his Ph.D. thesis [3]. Later on, he discovered [4] two criteria for 
detecting some useless reductions that made the Grobner bases method a practical tool to solve 
a wide class of problems in polynomial ideal theory and in many other research areas of science 
and engineering [5]. We refer to the monograph [2] for details on the theory of Grobner bases. 

The concept of comprehensive Grobner bases can be considered as an extension of these bases 
for polynomials over fields to polynomials with parametric coefficients. This extension plays an 
important role in application to constructive algebraic geometry, robotics, electrical network, au- 
tomatic theorem proving and so on (see, for example, [11-14]). Comprehensive Grobner bases and 
equivalent to them comprehensive Grobner systems were introduced in 1992 by Weispfenning [22]. 
He proved that any parametric polynomial ideal has a comprehensive Grobner basis and described 
an algorithm to compute it. In 2002, Montes [13] proposed a more efficient algorithm (DisPGB) 
for computing comprehensive Grobner systems. A year later Weispfenning in [21] proved the ex- 
istence of a canonical comprehensive Grobner basis. In 2003, Sato and Suzuki [17] introduced 
the concept of alternative comprehensive Grobner bases. Then in 2006, Manubens and Montes 
in [11] by using discriminant ideal improved DisPGB, and in [12] they introduced an algorithm 
for computing minimal canonical Grobner systems. Also in 2006, Sato and Suzuki [18] (see also 
[19]) suggested an important computational improvement for comprehensive Grobner bases by 
constructing the reduced Grobner bases in polynomial rings over ground fields. In 2010, Kapur, 
Sun and Wang [10], by combining Weispfenning's algorithm [22] with Suzuki and Sato's algorithm 
[18], proposed a new algorithm for computing comprehensive Grobner systems. More recently, in 
2010, Montes and Wibmer in [15] presented the GrobnerCover algorithm (its implementation in 
Singular is available at http://www-ma2.upc.edu/^montes/) which computes a finite partition 
of the parameter space into locally closed subsets together with polynomial data and such that 
the reduced Grobner basis for given values of parameters can immediately be determined from the 
partition. 

Involutive bases form an important class of Grobner bases. The theory of involutive bases goes 
back to the seminal works of French mathematician Janet. In the 20s of the last century, he de- 
veloped [9] a constructive approach to analysis of certain systems of partial differential equations 
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based on their completion to involution (cf. [20]). Inspired by the involution methods described 
in the book by Pommaret [16], Zharkov and Blinkov [23] introduced the concept of involutive 
polynomial bases in commutative algebra in the full analogy with the concept of involutive sys- 
tems of homogeneous linear partial differential equations with constant coefficients and in one 
dependent variable. Besides, Zharkov and Blinkov designed the first algorithm for construction 
of involutive polynomial bases. The particular form of an involutive basis they used is nowadays 
called Pommaret basis [20]. 

Gerdt and Blinkov [7] proposed a more general concept of involutive bases for polynomial 
ideals and designed efficient algorithmic methods to construct such bases. The underlying idea of 
the involutive approach is to translate the methods originating from Janet's approach into the 
polynomial ideals theory in order to provide a method for construction of involutive bases by 
combining algorithmic ideas in the theory of Grobner bases with constructive ideas in the theory 
of involutive differential systems. In doing so, Gerdt and Blinkov [7] introduced the concept of 
involutive division. Moreover, they derived the involutive form of Buchberger's criteria. This led 
to a strong computational tool which is a serious alternative to the conventional Buchberger 
algorithm. We refer to Seiler's book [20] for a comprehensive study and application of involution 
to commutative algebra and geometric theory of partial differential equations. 

In this paper, we introduce a notion of comprehensive involutive systems. For a parametric 
ideal, we decompose the space of parameters into a finite set of cells, and for each cell we yield the 
corresponding involutive basis of the ideal. Thereby, for each values of parameters, we find first a 
cell containing these values. Then, by substituting these values into the corresponding basis, we 
get the involutive basis of the given ideal. Based on the Gerdt-Blinkov involutive (abbreviated 
below by GBI) algorithm as described in [6] and also the Montes DisPGB algorithm [13], we 
present an algorithm for constructing comprehensive involutive systems. The proposed algorithm 
has been implemented in Maple, and we provide an illustrative example showing the step-by-stcp 
results of the algorithm. 

The paper is structured as follows. Section 2 contains the basic definitions and notations related 
to comprehensive Grobner systems, and a short description of the DisPGB algorithm. The basic 
definitions and notations from the theory of involutive bases are given in Section 3. In Section 4, 
the notion of comprehensive involutive system is introduced, and an algorithm for construction of 
such systems is described. In Section 5, we give an example illustrating in detail the performance 
of the algorithm of Section 4. 

2 Comprehensive Grobner Systems 

In this section, we recall the basic definitions and notations in theory of comprehensive Grobner 
systems and briefly describe the DisPGB algorithm. 

Let R = K[x] be a polynomial ring, where x = xi, . . . , x n is a sequence of variables and K is an 
arbitrary field. Below, we denote a monomial x" 1 ■ ■ ■ x" n <G R by x Q where a — (a\, . . . , a n ) S N™ 
is a sequence of non-negative integers. We shall use the notations deg^x") := Qj, deg(x a ) := 
Y^i—i a i- An admissible monomial ordering on R is a total order -< on the set of all monomials 
such that for any a, ft, 7 <G N n the following holds: 



A typical example of admissible monomial ordering is the lexicographical ordering, denoted 
by -<!iex- If x",x^ £ R arc two monomials, then x" -<!i ex x* 3 if the leftmost nonzero entry of 
ft — a is positive. Another typical example is the degree-reverse-lexicographical ordering denoted 
by ^dogroviex and defined as x a ^dogroviox x* 3 if dcg(x Q ) > deg(x /3 ) or dcg(x") = deg(x /3 ) and the 
rightmost nonzero entry of ft — a is negative. 

We shall write / = • ■ • , /fe) for the ideal I in R generated by the polynomials fi, ■ ■ ■ , fk € 
R. Let / e R and ^ be a monomial ordering on R. The leading monomial of / is the largest 
monomial (with respect to -<) occurring in /, and we denote it by LM(/). If F C R is a set of 
polynomials, then we denote by LM(F) the set {LM(/) | / € F} of its leading monomials. The 
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leading coefficient of /, denoted by LC(/), is the coefficient of LM(/). The leading term of / is 
LT(/) = LC(/)LM(/). The leading term ideal of I is defined as LT(J) = (LT(/) | / G I). 
A finite set G = {gi, . . . ,gk} C I is called a Grobner basis of I if LT(J) = (LT(gi), . . . , LT(g fc )}. 
For more details and definitions related to Grobner bases we refer to [2]. 

Now consider F = {/i, . . . , /&} C S := K[a, x] where a = ai, . . . , a m is a sequence of parame- 
ters. Let -< x (resp. -< a ) be a monomial ordering for the power products of Xi's (resp. <Zi's). We also 
need a compatible elimination product ordering -< x . a . It is defined as follows: For all a, 7 £ Z" 
and 0, 5 € Z^ 

x 7 a 5 -< x , a x. a a f3 x 7 -<; x x a or x 7 = x Q and a s ^ a a 13 . 

Now, we recall the definition of a comprehensive Grobner system for a parametric ideal. 

Definition 1. ([22]) A triple set {(Gi,Ni, Wj)}f =1 is called a comprehensive Grobner system for 
(F) w.r.t ^ X;a if for any i and any homomorphism a : K[a] — > K' (where K' is a field extension 
of K) satisfying 

(1) {VpGNiC K[a] ) [a(p) = 0], (ii) (Vq e W l C K[a] ) [a(q) ^ 0] 

we have cr(Gi) is a Grobner basis for o~((F)) C -^'[x] w.r.t. -< x . 

For simplification, we shall use the abbreviation CGS to refer to a comprehensive Grobner 
system, and CGSs in the plural case. For each i, the set Ni (resp. Wi) is called a (resp. non-) null 
conditions set. The pair (Ni, Wi) is called a specification of the homomorphism a if both conditions 
in the above definition are satisfied. 

Example 1. Let F — {ax 2 y — y 3 ,bx + y 2 } C K[a, 6, x, y] where a — a, b and x = x, y. Let us 
consider the lexicographical monomial ordering b ^i cx a on the parameters and on the variables 
y ^lox x as well. Using the DisPGB algorithm we can compute a CGS for (F) which is equal to 

{-&V + ay 5 ,bx + y 2 } {} {a, 6} 

{x 2 y,y 2 } {b} {a} 

{y 3 ,bx + y 2 } {a} {b} 

{y 2 } {a,b} {}. 

For instance, if a = 0, b = 2, then the third element of this system corresponds to this special- 
ization. Therefore, {y 3 , 2x + y 2 } is a Grobner basis for the ideal (F)| a =o.b=2 = (— J/ 3 , 2a; + y 2 ). 

Remark that, by the above definition, a CGS is not unique for a given parametric ideal, and one 
can find other partitions for the space of parameters, and, therefore, other CGSs for the parametric 
ideal. 

Now, we briefly describe the Montes DisPGB algorithm to compute CGSs for parametric ideals 
(see [13, 11]). The main idea of DisPGB is based on discussing the nullity or not w.r.t. a given 
specification (N, W) for the leading coefficients of the polynomials appearing at each step (this 
process is performed by the NewCond subalgorithm). Let us consider a set F C £ of parametric 
polynomials. Given a polynomial / e F and a specification (N,W), NewCond is called. Three 
cases are possible: If LC(/) specializes to zero w.r.t. (N, W), we replace / by / — LT(/), and then 
start again. If LC(/) specializes to a nonzero element we continue with the next polynomial in F . 
Otherwise (if LC(/) is not decidable, i.e. we can't decide whether or not it is null w.r.t. (iV, W)), 
the subalgorithm Branch is called to create two complementary cases by assuming LC(/) = 
and LC(/) 7^ 0. Therefore, two new disjoint branches with the specifications (N U {LC(/)}, W) 
and (N, W U {LC(/)}) are made. This procedure is continued until every polynomial in F has a 
nonnull leading coefficient w.r.t. the current specification. Then, we proceed with CondPGB: This 
algorithm receives, as an input, a set of parametric polynomials and a specification (N, W) and, by 
applying Buchbcrger's algorithm, creates new polynomials. When a new polynomial is generated, 
NewCond verifies whether its leading coefficient leads to a new condition or not. If a new condition 
is found, then the subalgorithm stops, and Branch is called to make two new disjoint branches. 
Otherwise, the process is continued and computes a Grobner basis for (F) , according to the current 
specification. The collection of these bases, together with the corresponding specifications yields 
a CGS for (F). 
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3 Involutive Bases 



Now we recall the basic definitions and notations concerning involutive bases and present below 
the general definition of involutive bases. First of all, we describe the cornerstone notion of invo- 
lutive division [7] as a restricted monomial division [6] which, together with a monomial ordering, 
determines properties of an involutive basis. This makes the main difference between involutive 
bases and Grobner bases. The idea behind involutive division is to partition the variables into two 
subsets of multiplicative and nonmultiplicative variables, and only the multiplicative variables can 
be used in the divisibility relation. 

Definition 2. [7,6] An involutive division C on the set of monomials of R is given, if for any 
finite set U of monomials and any u G U , the set of variables is partitioned into subsets Mc{u, U) 
of multiplicative and N Mc(u,U) of nonmultiplicative variables such that 

1. u, v G U, uC{u, U) n vjC(v, U) ^ => u G vC(v, U) or v G uC(u, U), 

2. v G U, v G uC(u, U) => C(v, U) C C(u, U), 

3. ueV andV CU=^ C{u, U) C C(u, V), 

where C(u, U) denotes the set of all monomials in the variables in Mc(u, U).IfvE uC(u, U), then 
we call u an C— (involutive) divisor of v, and we write u\cv. If v has no involutive divisor in a set 
U, then it is £— irreducible modulo U . 

In this paper, we are concerned with the wide class [8] of involutive divisions determined by 
a permutation p on the indices of variables and by a total monomial ordering □ which is either 
admissible or the inverse of an admissible ordering. This class is defined by 

(Vuetf) [NM^(u,U)= (J NM^(u,{u,v}) ] (1) 

veu\{u} 

where 

{if u □ v or (u □ v A v I u) then 
else {x pW }, i = minjj | deg pU) (u) < deg p(j) {v)\ . 

Remark 1. The involutive Janet division introduced and studied in [7] is generated by formulae 
(l)-(2) if □ is the lexicographic monomial ordering ^i cx and p is the identical permutation. The 
partition of variables used by Janet himself [9] (see also [20]) is generated by >-i ex as well with the 
permutation which is inverse to the identical one: 

I 2 ... n 
n n — 1 ... 1 

Throughout this paper £ is assumed to be a division of the class (l)-(2). Now, we define an 
involutive basis. 

Definition 3. Let I C R be an ideal, -< be a monomial ordering on R and C be an involutive 
division. A finite set G C I is an involutive basis of I if for all f G J there exists g G G such that 
LM(g)\cLM(f). An involutive basis G is minimal if for any other involutive basis G the inclusion 
LM(G) C LM(G) holds. 

^From this definition and from that for Grobner basis [3, 2] it follows that an involutive basis 
of an ideal is its Grobner basis, but the converse is not always true. 

Remark 2. By using an involutive division in the division algorithm for polynomial rings, we obtain 
an involutive division algorithm. If G is an involutive basis for an involutive division C, we use 
NF/;(/, G) to denote C— normal form of / modulo G, i.e. the remainder of / on the involutive 
division by G. A polynomial set F is C— autoreduced if / = NF,c(/, F \ {/}) for every / G F. 
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The following theorem provides an algorithmic characterization of involutive bases which is an 
involutive analogue of the Buchberger characterization of Grobner bases. 

Theorem 1. ([7,8]) Given an ideal I C R, an admissible monomial ordering -< on R and an 
involutive division C, a finite subset G d I is an involutive basis of I if for each f G G and each 
x G N Mc(LM(f) ,LM(G)) the equality NFc( x f,G) = holds. An involutive basis exists for any 
I, C and -<. A monic and C-autoreduced involutive basis is uniquely defined by I and -<. 

4 Comprehensive Involutive Systems 

In this section, like the concept of comprehensive Grobner systems, we define the new notion of 
comprehensive involutive system for a parametric ideal. Then, based on the GBI algorithm [6] 
and the Montes DisPGB algorithm [13], we propose an algorithm for computing comprehensive 
involutive systems. 

Definition 4. Consider a finite set of parametric polynomials F C S — K[&,x] where K is a 
field, x = xi, . . . , x n is a sequence of variables and a = <Zi, . . . , a rn is a sequence of parameters, -< x 
(resp. -< a ) is a monomial ordering involving the Xi's (resp. ai's), and C is an involutive division 
on K\x\. Let M — {(Gi, Ni, Wj)}f =1 be a finite triple set where sets Ni,Wi C K[a] and G, C S 
are finite. The set M is called an (C— ) comprehensive involutive system for (F) w.r.t -<; x ,a if for 
each i and for each homomorphism a : K[a] — > K 1 (where K' is a field extension of K) satisfying 

(i) (VpeJVi) [<Kp)=0], (ii) (VqeWi) [o-(q)^0] 

a(Gi) is an (C-)involutive basis for a ((F)) c K'\x\. We use the abbreviation CIS (resp. CISs) to 
stand for comprehensive involutive system (resp. systems). M is called minimal, if for each i, the 
set o~(Gi) is a minimal involutive basis. 

Given a CGS, one can straightforwardly compute a CIS by using the following proposition. 

Proposition 1. Let G = {gi, . . . , gt} be a minimal Grobner basis of an ideal I c K[x\, . . . , x n ] 
for a monomial ordering -<. Let hi = maxggcKdeg^LMG?))}. Then the set of products 

{mg | g e G, m is a monomial s.t. (Vi) [degj(ra) < hi — degj(LM(g)) ]} (3) 

is an C-involutive basis of I . 

Proof. Denote LM(G) by U. From (1) (2) it follows 

(VueU) (Vx t e NM c (u,U))) [Acg t (u)<K]. (4) 

It is also clear that if we enlarge G with a (not necessarily nonmultiplicative) prolongation gxj of 
its element g G G such that deg J (LM(g)) < hj, then (4) holds for the enlarged leading monomial 
set U := U U {LM(g)xj} as well. Consider completion G of the polynomial set G with all possible 
prolongations of its elements satisfying (3) and denote the monomial set LM(G) by U. Then 

(Vu G U) (Vx G NM c (u,U)) (3v G U) [v \ c ux] . 

This means, by Theorem 1, that the monomial set U is an involutive basis of (LM(G)). Now, since 
G is a Grobner basis of I we have LT(7) = (LM(G)), and hence LT(7) = (LM(G)>. Therefore, G 
is an involutive basis of I by Definition 3. □. 

Example 2. Let F — {ax 2 , by 2 } G K[a, x] where a — a, b and x = x, y. Let also b <i ex a and 
V <iex x. Then, F is a CGS for any sets of null and nonnull conditions. Using Proposition 1, we 
can construct the following Janet basis of (F) which is a GIS for any sets of null and nonnull 
conditions: 

r 2i2 2 22?. 2j2 2i 

{ax , by , ayx ,ay x , bxy ,bxy\. 
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On the other hand, the algorithm that we present below computes the following minimal Janet 
CIS for (F): 

{ax 2 ,by 2 ,bxy 2 } {} {a, b} 
{ax 2 } {b} {a} 

{by 2 } {a} {b} 

{0} {a, 6} {}. 

Remark 3. Using Proposition 1, we cannot directly compute a minimal CIS from a given CCS. 
Indeed, to do this, we must examine the leading coefficients of each Grobner basis in the CGS, and 
this may lead to further partitions of the space of parameters. Moreover, the CIS computed by this 
way may be too large, since many prolongations constructed by means of (3) may be useless. That 
is why, based on the GBI algorithm [6] and on the Montes DisPGB algorithm [13], we propose a 
more efficient algorithm for computing minimal CISs. 

Now we describe the structure of polynomials that is used in our new algorithm. To avoid 
unnecessary reductions (during the computation of involutive bases) by applying the involutive 
form of Buchberger's criteria (see [6]), we need to supply polynomials with additional structural 
information. 

Definition 5. [6] An ancestor of a polynomial f £ F C i?\{0}, denoted by anc(/) 7 is a polynomial 
g E F of the smallest deg(LM(g)) among those satisfying LM(/) = uLM(g) where u is either 
the unit monomial or a power product of nonmultiplicative variables for LM(g) and such that 
NF c (f-ug,F\{f}) = Oiff^ug. 

Below we show how to use this concept to apply the involutive form of Buchberger's criteria. In 
what follows, we store each polynomial / as the p = [/, g 7 V] where / = poly(p) is the polynomial 
part of p, g = anc(p) is the ancestor of / and V = NM(p) is the list of nonmultiplicative 
variables of / have been already used to construct prolongations of / (sec the for-loop 20-23 in 
the subalgorithm GBI). If P is a set of triples, we denote by poly(P) the set {poly(p) | p E P}. 
If no confusion arises, we may refer to a triple p instead of poly(p), and vice versa. 

We present now the main algorithm ComInvSys which computes a minimal CIS for a given 
ideal. It should be noted that we use the subalgorithms NewCond and CanSpec (resp. Tail- 
NormalForm) as they have (resp. it has) been presented in [13] (resp. [6]), and recall them for 
the sake of completeness. Also, we use the subalgorithm Branch (resp. GBI , HeadReduce and 
HeadNormalForm) from [13] (resp. [6]) with some appropriate modifications. 

Algorithm ComInvSys 



Input: F, a set of polynomials; £, an involutive division; -< x , a monomial ordering on the 

variables; -< a , a monomial ordering on the parameters 
Output: a minimal CIS for (F) 

1: global: List, ind; 

2: List:=Null; 

3: ind:=l; 

4: B:={[F[t\,F[i\,HS\ M = 1,...,|F|}; 

5: G:={Branch([F[1],F[1],0],B,{ },{ },{ })}; 

6: for i from 2 to |F| do 

7: ind:=ind+l; 

8: G:={BRANCu([F\i},F[i\,<b],A[2},A[3},A[A],A{5]) A E G}; 
9: od 

10: Return (List) 



In the above algorithm, List is a global variable to which we add any computed involutive basis 
together with its corresponding specification to form the final CIS. That is why, at the beginning 
of computation we must set it to the empty list (see Branch). Note that here and in Branch, we 
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use |F| to denote the number of polynomials in the input set F. The variable ind is also a global 
variable, and we use it to examine all the leading coefficients of the elements in F (see Branch). 
Once we are sure about the non-nullity of these coefficients, then we start the involutive basis 
computation. Indeed, Branch inputs a triple p = [f,g,V], a set B of examined and processed 
polynomials, a set N of null conditions, a set W of nonnull conditions and a set P of non-processed 
polynomials. Then, it analyses the leading coefficient of / w.r.t. N and W. Now, two cases are 
possible: 

— ind< \F\: If LC(/) is not decidable by N and W then we create two complementary cases by 
assuming LC(/) = and LC(/) ^ 0. Then we pass to the next polynomial in F. 

— ind= \F\: We are now sure that we have examined all the leading coefficients of the elements 
in F (except possibly the very last one which is to be /). If LC(/) is not decidable by N and 
W then we again create two complementary cases with LC(/) = and LC(/) ^ 0. Otherwise, 
we continue to process the polynomials in P by using the GBI algorithm. If P = this means 
that B is an involutive basis consistent with the conditions in N and W, and we add (B, N, W) 
to List. 

It is worth noting that if the input specification of Branch is incompatible, then it stops 
the process only for the corresponding branch, and continues the construction of other branches. 
Moreover, using the above notations, if ind< |F| and no new condition is detected, then Branch 

returns an element of the form (p,B N ,N',W',P) where p is a triple, N',W are two sets of 

— N' 

conditions, B is the normal form of a specializing basis B and P is a set of non-examined triples. 
Otherwise, it calls itself to create the new branches. Finally, if ind= \F\, then the algorithm does 
not return anything and completes the global variable List. 

Subalgorithm Branch 

Input: p, a triple; B, a specializing basis; N, a set of null conditions; W, a set of nonnull 

conditions; P, a set of non-examined triples 
Output: It stores the refined (B' , N', W , P'), and creates two new vertices when necessary 

or marks the vertex as terminal 
1: P=[f,9,V\; 

2: (test, N,W):=CanSpec(N,W); 
3: if test=false then 

4: Return STOP (incompatible specification has been detected) 
5: fi 

6: (cd,f',N',W) :=NewCond(/, iV, w): 

7: p := [f',g N , V] (g N denotes the remainder of the division of g by N'); 
8: if ind < \F\ and cd ^ { } then 

9: BRANCu(p,B,N',W'l)cd,P); Branch (p, B, N' U cd, W, P); 
10: fi 

11: if ind < \F\ and cd = { } then 

12: Return (p,B~ N ,N',W',P) 
13: fi 

14: if cd = { } then 

15: {test,p', B', N', W, P') :=GBI (B, N', W, P); 

16: if test then 

17: List .-List, (B',N',W); 

18: else 

19: BRANCH(p',B',iV',VK',P'); 
20: fi 
21: else 

22: BRANCU(p, B,N',W U cd, P); Branch(p, B, N' U cd, W, P); 
23: fi 
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The subalgorithm CanSpec produces a quasi-canonical representation for a given specification. Its 
subalgorithm FacVar invoked in lines I and 13 returns the set of factors of its input polynomial. 

Definition 6. ([13]) A specification (N,W) is called quasi- canonical if 

— TV is the reduced Grobner basis w.r.t. -< a of the ideal containing all polynomials that specialize 
to zero in K [a] . 

— The polynom ials in W specializing to non-zero are reduced modulo TV and irreducible over K[&] 



Input: TV, a set of null conditions; W, a set of nonnull conditions 

Output: true if TV and W are compatible and false otherwise; (TV', W), a quasi-canonical 
representation of (TV, W) 
1: W :=FACVAR({g Ar : q G W}); test -true; TV' := TV; h := ]J qeW q; 

2: if h G ^f{N J ) then 
3: iesi:=false; TV' := {1}; 
4: Return (test, N' ,W); 
5: fl 

6: flag:— true; 
7: while flag do 
8: flag:=false; 

9: N":= Remove any factor of a polynomial in N' that belongs to W; 
10: if N" / N' then 
11: flag:=tine; 

12: N':= a Grobner basis of (N") w.r.t. ^ a ; 

13: W :=FACVAR({g Ar ' : q G W'})\ 

14: fi 
15: od 

16: Return (test, N' , W) 



Input: /, a parametric polynomial; TV, a set of null conditions; W, a set of nonnull condi- 
tions 

Output: cd, a new condition; /', a parametric polynomial; N', a set of null conditions; W, 
a set of nonnull conditions 
1: /' := /; test:=true; N' := N; cd:={ }; 
2: while test do 
3: if LC(/') G yW) then 

4: TV' := a Grobner basis for (TV', LC(/')) w.r.t. ^ a ; 
5: /' := /' - LT(/); 

6: else 

7: test:=f&\se; 
8: fi 
9: od 

10: /' := /' ; 

11: W := {W N ' | w G W}; 

12: cd := cd U FacVar(LC(/')) \ W; 

13: Return(cd, /', TV', W) 




Subalgorithm CanSpec 



Subalgorithm NewCond 



Comprehensive Involutive Systems 9 



We describe now the NewCond subalgorithm. When it is invoked in line 6 of Branch with the 
input data (/, N, W), one of the two following cases may occur: 

1. If LC(/) is dccidablc w.r.t. the specification (N, W), then the subalgorithm returns: 

(i) NewCond(/ - LT(f),N, W) in the case when LC(/) specializes to zero w.r.t. (N, W). 

(ii) (0, /, N, W) in the case when LC(/) does not specialize to zero w.r.t. (N, W). 

2. If LC(/) is not decidable w.r.t (N,W), then NewCond returns (cd,f,N,W) where set cd 
contains one of the non-decidable factors (w.r.t (N, W)) of LC(/). 

It should be emphasized that FacVar(LC(/')) \ W in line 12 returns only one factor of LC(/'). 

The subalgorithm GBI, presented below, is an extension of the algorithm InvolutiveBasis 
II described in [6]. The latter algorithm computes involutive bases and applies the involutive form 
of Buchbcrgcr's criteria to avoid some unnecessary reductions [7] (see also [1,6]). The criteria are 
applied in the subalgorithm HeadNormalForm (see line 7) that is invoked in line 5 of GBI. 

Proposition 2. ([7, 6]) Let I C R be an ideal and G G I be a finite set. Let also -< be a monomial 
ordering on R and C be an involutive division. Then G is an £— involutive basis of I if for all 
f G G and for all x G NMc(LM(f ), LM(G)) one of the two conditions holds: 

1. m c {xf,G) = o. 

2. There exists g G G with LM(g)|£LM(x/) satisfying one of the following conditions: 
(Ci) LM(anc(/))LM(anc( 5 )) = LM(xf) , 

(C 2 ) lcm(LM(anc(/)), LM(anc(g))) is a proper divisor ofLM(xf) . 
Subalgorithm GBI 

Input: B, a specializing basis; N, a set of null conditions; W, set of nonnull conditions; P, 

set of non-examined triples 
Output: If tesi=true, a minimal involutive basis for (B) w.r.t. C and -< x ,a; otherwise, it 

returns a triple so that we must discuss the leading coefficient of its polynomial part 

1: if P = then 

2: Select p G B with no proper divisor of LM(poly(p)) in LM(poly(B)) 

3: T:={p}; Q:=B\{p}; 

4: else 

5: V : l>>: Q := P; 

6: fi 

7: while Q ^ do 

8: {test,p,T,N,W,Q') :=HeadReduce(T, N, W, Q)\ 

9: if test =false then 

10: Return (false, p, T, N, W, Q') 

11: fi 

12: Q := Q'\ 

13: Select and remove p G Q with no proper divisor of LM(poly(p)) in LM(poly(Q)); 

14: if poly(p) = anc(p) then 

15: for q G T whose LM(poly(g r )) is a proper multiple of LM(poly(p)) do 

16: Q:=QU{q}; T:=T\{q}; 

17: od 

18: fi 

19: h :=TailNormalForm(/>, T); T := T U {{h, anc(p), NM(p)}}; 

20: for q G T and x G ^M £ (LM(poly(q)), LM(poly(T))) \ NM(q) do 

21: Q:=Qu{{xpoly(9),anc( g ),0}}; 

22: NM(q) := NM(q) fl JVM £ (LM(poly(g)), LM(poly(T))) U {x}; 

23: od 

24: od 

25: Return (true, 0, T, N, W, { }) 
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This algorithm invokes three subalgorithms HeadReduce, TailNormalForm and HeadNor- 
MALForm that we present below. The subalgorithm HeadReduce performs the involutive head 
reduction of polynomials in the input set of triples modulo the input specializing basis. The subal- 
gorithm TailNormalForm (resp. HeadNormalForm) invoked in line 19 of GBI (rcsp. in line 
4 of HeadReduce) computes the involutive tail normal form (resp. the involutive head normal 
form) of the polynomial in the input triple modulo the input specializing basis. 

In the subalgorithm HeadNormalForm, the Boolean expression Criteria(p, g) is true if at 
leat one of the conditions (Ci) or (C2) in Proposition 2 are satisfied for p and g, false otherwise. 
We refer to [6] for more details on the algorithm GBI and on its subalgorithms. 



Subalgorithm HeadReduce 



Input: B, a specializing basis; N, a set of null conditions; W, a set of nonnull conditions; 




P a set of non-examined triples 


Output: If test=true, the £-head reduced form of P modulo B; otherwise, it returns a 




triple such that we must examine the leading coefficient of its polynomial part 


1 


S:=P; g:=0; 


2 


while S ^ do 


3 


Select and remove peS; 


4 


(test, h, B, N, W) :=HeadNormalForm(p, B, N, W); 


5 


if test=ialse then 


6 


Return (false, p, B, N,W,SU Q) 


7 


fi 


8 


if h ^ then 


9 


if LM(poly(p)) ^ LM(ft) then 


10 


Q :=QU{{h,h,9}}; 


11 


else 


12 


Q:=QU{p}; 


13 


fi 


14 


else 


15 


if LM(poly(p)) = LM(anc(p)) then 


16 


S:=S\{qeS | anc(g) = poly(p)}; 


17 


fi 


18 


fi 


19 


od 


20 


Return (true, 0, B, N, W, Q) 



Subalgorithm TailNormalForm 



Input: p, a triple; B, a set of triples 


Output: ^-normal form of poly(p) modulo poly(B) 


1 


h := poly(p); 


2 


G := poly(B); 


3 


while h has a term t which is £— reducible modulo G do 


4 


Select g e G with LM(g)\ c t: 


5 


h .— h LT(g) ' 


6 


od 


7 


Return (h) 
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Subalgorithm HeadNormalForm 

Input: p, a triple; B, a specializing basis; N, a set of null conditions; W, set of nonnull 
conditions 

Output: If test=truc, the £-head normal form of poly(p) modulo B; otherwise, a polyno- 
mial whose leading coefficient must be examined 
1: ft := poly(p); G := poly(B); 
2: if LM(/i) is ^-irreducible modulo G then 
3: Return (true, ft, B, N, W) 
4: else 

5: Select g G G with LM(poly(ff))| £ LM(/i); 
6: if LM(h) ^ LM(anc(p)) then 
7: if Criteria(p, g) then 
8: Return (true, 0, B, N, W) 

9: fi 
10: else 

11: while ft 7^ and LM(ft) is /^-reducible modulo G do 
12: Select 5 e G with LM(3)| £ LM(/i); 

13: ft:=ft-<?gjg}; 

14: (erf, ft', AT', W') :=NEWCOND(/l, AT, W); 

15: if erf ^ then 

16: Return (false, ft', B, AT', IF') 

17: fi 

18: od 

19: fi 

20: fi 

21: Return (true, ft, £?, AT, W) 



Theorem 2. Algorithm ComInvSys terminates in finitely many steps, and computes a minimal 
CIS for its input ideal. 

Proof. Let I — (F) where F = {fi, ■ ■ ■ , fk} C AT[a,x] is a parametric set, x = x\,...,x n (resp. 
a = ai, . . . , a m ) is a sequence of variables (resp. parameters). Let -< x (resp. -< a ) be a monomial 
ordering involving the Xi's (resp. a^'s), and C be an involutive division on AT[x]. 

Suppose that ComInvSys receives F as an input. To prove the termination, we use the fact 
that AT [a] is a Noetherian ring. When Branch is called, the leading coefficient of some polynomial 
/ G I is analyzed. For this purpose, the subalgorithm NewCond determines whether LC(/) is 
decidable or not w.r.t. the given specification (N,W). Two alternative cases can take place: 

— LC(/) is decidable and we check the global variable ind. Now if ind< k, then we study the 
next polynomial in F. Otherwise, GBI is called. If all the leading coefficients of the examined 
polynomials (to compute a minimal involutive basis) are decidable, then the output, say G, is 
a minimal involutive basis of I w.r.t. (AT, W), and we add (G, N, W) to List. Otherwise, two 
new branches are created by calling Branch (cf. the second case given below). In doing so, 
the minimality of G and the termination of its computation is provided by the structure of 
GBI algorithm (see [6]). 

— LC(/) is not decidable and we create two branches with (N, W U erf) and (N U erf, W), where 
erf is the one-element set containing the new condition derived from LC(/). 

Thus, in the second case, the branch for which A^ (resp. W) is assumed, increases the ideal 
(N) C AT[a] (resp. (W) C AT[a]). Note that we replace N by a Grobner basis of its ideal (see line 
4 in NewCond). Since the ascending chains of ideals stabilize, the algorithm terminates. This 
argument was inspired by the proof in [13], Theorem 16. 

To prove the correctness, assume that M — {(Gj, AT*, Wj)}f =1 is the output of ComInvSys 
for the input is F (note that we have used the fact the this algorithm terminates in finitely many 
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steps). Consider integer 1 < i < I homomorphism u : K[a] — > K' where (Ni, W{) is a specification 
of o and K' is a field extension of K. 

We have to show that for each fed and each x E NMc(LM(a(f)),LM(a(Gi))), in accor- 
dance with Theorem 1, the equality NFc(cr(xf), (7(d)) = holds. By using 'reductio ad absur- 
dum', suppose g — NF£(ct(x/), <j{Gi)) and g ^ 0. Since (Gi, Ni, Wj) has been added to List in 
Branch, the leading coefficients of the polynomials in the subalgorithm GBI, examined at com- 
putation of a minimal involutive basis for F, are decidable w.r.t. (Ni, Wi). Furthermore, / 6 Gi 
implies that in the course of GBI xf is added to Q, the set of all nonmultiplicative prolongations 
that must be examined (see the notations used in GBI). Then, HeadReduce is called to perform 
the £-head reduction of the elements of Q modulo the last computed basis T C Gi. The computed 
£-head normal form of xf is further reduced by invoking TailNormalForm which performs the 
£-tail reduction. By the above notations, g is the result of this step. Thus, g should be added to 
T C Gi. It follows that NFc(a(xf),a(Gi)) — 0, a contradiction, and this completes the proof. □ 

5 Example 

Now we give an example to illustrate the step by step construction of a minimal CIS by the 
algorithm ComInvSys proposed and described in the previous section 3 . 

For the input F = {ax 2 , by 2 } C K[a,b,x,y] from Example 2, Janet division and the lexico- 
graphic monomial ordering with b -<\ cx a and y -<\ eK x the algorithm performs as follows: 

^COMlNVSYSfT, C, ^lex, <lex) 

List := Null; ind := 1; k := 2; 

B := {[ax 2 , ax 2 , 0], [by 2 , by 2 ,0]} 
^BRANCH([ax 2 ,ax 2 ,0],S,{ },{ },{ }) 

^NEwCOND(ax 2 ,{ },{}) = ({a},{ },{ }) 

^BRANCH([ax 2 ,ax 2 ,0],S,{ },{a},{ }) 

->NEwCOND(ax 2 , { }, {a}) = ({},{ }, {a}) 
G := {{{ax 2 ,ax 2 ,%B,{},{a},{})} 

^BRANCH([ax 2 ,ax 2 ,0],S,{a},{ },{ }) 

^NEwCOND(ax 2 ,{ a },{ }) = ({ },{a},{ }) 
G := { ([ax 2 , ax 2 , 0], B, { }, {a}, { }), ([ax 2 , ax 2 , 0], {[0, 0, 0], [by 2 , by 2 , 0]}, {a}, {},{})} 
ind := 2; 

A = ([ax 2 , ax 2 , %B, {},{«},{}) 

— >Branch([&2/ 2 , by 2 , 0], _B, { },{a},{ }) 

^NEwCOND(6y 2 , { }, {a}) = ({b}, {},{}) 
^BRANCH([6 J / 2 ,6y 2 ,0], J B,{ },{a,b},{ }) 

(* further Br AN en ([by 2 , by 2 , ®},B,{b}, {a}, { }) is executed*) 
-^NEwCOND( y 2 ,{ },{a,b}) = ({ },{ },{a,b}) 
ind > k = 2 
cd = {} 

^GBI (B,{ },{a,b},{ }) 
T:={[by 2 ,by 2 ,m 
Q :={[ax 2 ,ax 2 ,0]} 
^HeadReduce(T, { }, {a, b}, Q) 

->HEADNoRMALFORM([ax 2 ,ax 2 ,0],T, { },{a,b}) = (true, ax 2 ,T, { },{a, 6}) 
HeadReduce returns (true, 0, T, { }, {a, b}, Q) 
p := [ax 2 ,ax 2 ,0] 

3 The Maple code of our implementation of the algorithm for the Janet division [7] is available at 
http : // invo . j inr . ru and http : / / amirhashemi . iut . ac . ir/ software . html 
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Q = {} 

-^TailNormalForm(j>, T) = ax 2 
T ■= {[by 2 , by 2 , 0], [ax 2 , ax 2 , 0]} 
Q := {[bxy 2 ,by 2 ,%]} 

-^HeadReduce(T, { }, {a, b}, Q) = {true, 0, T, { }, {a, b}, Q) 
p := [bxy 2 , by 2 , 0] 
Q = {} 

-^TailNormalForm(p, T) = bxy 2 
T := {[by 2 , by 2 , 0], [ax 2 , ax 2 , 0], [bxy 2 , by 2 , 0]} 
Q:={[bx 2 y 2 ,by 2 M 

-^HeadReduce(T, { }, {a, b}, Q) = (true, 0, T, { }, {a, b}, { }) 

Q--={} 

->GBI returns (true, 0, {by 2 , ax 2 , bxy 2 }, { }, {a, b}) 
List := ({by 2 , ax 2 , bxy 2 }, { },{a,b}) 
B = { [ax 2 , ax 2 , 0], [0,0,0]} 

— >Branch([6j/ 2 , by 2 , 0], B, {b}, {a}, { }) 

^NewCond(^ 2 , {b}, {a}) = ({ }, {b}, {a}) 
ind > k = 2 
cd = { } 

^GBI (B, {b}, {a}, { }) = (true, 0, {ax 2 }, {b}, {a}) 
List := ({by 2 , ax 2 , bxy 2 }, { }, {a, b}), ({ax 2 }, {b}, {a}) 
(* Return back to ComInvSys *) 
A = ([ax 2 , ax 2 , 0], {[0, 0, 0], [by 2 , by 2 ,%]}, {a}, {},{}) 
B = {[Q,i),%[by 2 ,by 2 ,%]} 
^BKAKCn([by 2 ,by 2 ,%],B,{ },{a},{ }) 

^NewCond(&2/ 2 , {a}, { }) = ({&}, {},{}) 
— >Branch([6j/ 2 , by 2 , %B, {a}, {&}, { }) 

(* further Branch([6j/ 2 , by 2 , 0], B, {a, b}, { },{ }) is executed *) 
->NEwCOND(6y 2 , {a}, {6}) = ({ }, {a}, {6}) 
ind > k = 2 
ed = {} 

^GBI (B, {a}, {6}, { }) - (true, 0, {6y 2 }, {a}, {6}) 
List := ({by 2 , ax 2 , bxy 2 }, { }, {a, &}), ({ax 2 }, {&}, {a}), ({6y 2 }, {a}, {6}) 
£J = {[0,0,0], [0,0,0]} 

^BRANCH([6j/ 2 ,6 J/ 2 ,0], J B,{a,6},{ },{ }) 

-^NEwCOND(foy 2 ,{a,6},{ }) = ({ },{a,b},{ }) 
ind > k = 2 
cd = {} 

^GBI (B, {a, &}, {},{}) = (true, 0, {0}, {a, b}, { }) 
List := ({by 2 , ax 2 , bxy 2 }, { }, {a, b}), ({ax 2 }, {b}, {a}), ({by 2 }, {a}, {&}), ({0}, {a, b}, { }). 



Acknowledgements 

The main part of research presented in the paper was done during the stay of the second author 
(A.H.) at the Joint Institute for Nuclear Research in Dubna, Russia. He would like to thank the 
first author (V.G.) for the invitation, hospitality, and support. The contribution of the first author 
was partially supported by grants 01-01-00200, 12-07-00294 from the Russian Foundation for Basic 
Research and by grant 3802.2012.2 from the Ministry of Education and Science of the Russian 
Federation. 



14 Vladimir Gcrdt and Amir Hashemi 



References 

1. Apel, J., Hemmecke, R.: Detecting unnecessary reductions in an involutive basis computation. J. 
Symbolic Computation 40, 1131-1149 (2005) 

2. Becker, T., Weispfcnning, T.: Grobner Bases: a Computational Approach to Commutative Algebra. 
Graduate Texts in Mathematics, 141, Springer- Verlag, New York (1993) 

3. Buchberger, B.: Ein Algorithms zum Auffinden der Basiselemente des Restklassenrings nach einem 
nuildimensionalen Polynomideal. PhD thesis, Universitat Innsbruck (1965) 

4. Buchberger, B.: A Criterion for Detecting Unnecessary Reductions in the Construction of Grobner 
Bases. In: Edward W. Ng (ed.) EUROSAM'79. LNCS, vol.72, pp. 3-21. Springer, Berlin (1979) 

5. Buchberger, B., Winkler, F. (eds.): Grobner Bases and Applications. London Mathematical Society 
Lecture Note Series, Vol. 251. Cambridge University Press, Cambridge (1998) 

6. Gerdt, V.P.: Involutive Algorithms for Computing Grobner Bases. In: Cojocaru, S. Pfister, G., Uf- 
narovski, V. (eds.) Computational Commutative and Non-Commutative Algebraic Geometry, pp. 199- 
225. IOS Press, Amstrerdam (2005) (arXiv:math/0501111) 

7. Gerdt, V.P., Blinkov, Yu.A.: Involutive Bases of Polynomial Ideals. Mathematics and Computers in 
Simulation 45, 519-542 (1998) 

8. Gerdt, V.P., Blinkov, Yu.A.: Involutive Division Generated by an Antigraded Monomial Ordering. In: 
Gerdt, V.P., Koepf, W., Mayr, E., Vorozhtsov, E.V. (eds.) CASC 2011. LNCS, vol. 6885, pp. 158-174. 
Springer, Berlin (2011) 

9. Janet, M.: Les Systemes d'Equations aux Derivees Partielles. Journal de Mathematique 3, 65-151 
(1920) 

10. Kapur, D., Sun, Y., Wand, D.: A New Algorithm for Computing Comprehensive Grobner Systems. 
In: Watt, S.M. (ed.) Proc. ISSAC'10, pp. 29-36. ACM Press, New York (2010) 

11. Manubens, M., Montes, A.: Improving DisPGB algorithm using the discriminant ideal. J. Symbolic 
Computation 41, 1245-1263 (2006) 

12. Manubens, M., Montes, A.: Minimal Canonical Comprehensive Grobner Systems. J. Symbolic Com- 
putation 44, 463-478 (2009) 

13. Montes, A.: A new algorithm for discussing Grobner bases with parameters. J. Symbolic Computation 
33, 183-208 (2002) 

14. Montes, A.: Solving the load flow problem using Grobner bases. SIGSAM Bulletin 29, 1-13 (1995) 

15. Montes, A., Wibmer, M.: Grobner bases for polynomial systems with parameters. J. Symbolic Com- 
putation 45, 1391-1425 (2010) 

16. Pommaret, J.-F.: Systems of Partial Differential Equations and Lie Pseudogroups. Mathematics and 
its Applications, vol. 14. Gordon & Breach Science Publishers, New York (1978) 

17. Sato, Y., Suzuki, A.: An alternative approach to comprehensive Grobner bases. J. Symbolic Compu- 
tation 36, 649-667 (2003) 

18. Sato, Y., Suzuki, A.: A simple algorithm to compute comprehensive Grobner bases using Grobner 
bases. In: Trager, B.M. (ed.) Proc. ISSAC'2006, pp. 326-331. ACM Press, New York (2006) 

19. Suzuki, A.: Computation of full comprehensive Grobner bases. In: Ganzha, V.G, Mayr, E., Vorozhtsov, 
E.V. (eds.) CASC 2005. LNCS, vol. 3781, pp. 431-444. Springer, Berlin (2005) 

20. Seiler, W.M.: Involution - The Formal Theory of Differential Equations and its Applications in Com- 
puter Algebra. Algorithms and Computation in Mathematics, vol. 24. Springer- Verlag, Berlin (2010) 

21. Weispfcnning, V.: Cannonical comprehensive Grobner bases. J. Symbolic Computation 36, 669-683 
(2003) 

22. Weispfenning, V.: Comprehensive Grobner Bases. J. Symbolic Computation 14, 1-29 (1992) 

23. Zharkov, A.Yu., Blinkov, Yu.A.: Involutive approach to investigating polynomial systems. Mathemat- 
ics and Computers in Simulation 42, 323-332 (1996) 



